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Abstract 

In this paper, a study of light-cone distribution ampHtudes for p-wave heavy quarkonium states 
are presented. Within the hght-front framework, the leading twist light-cone distribution ampli- 
tudes, and their relevant decay constants, have some simple relations. These relations can be 
further simplified when the non-relativistic limit and the wave function as a function of relative 
momentum |k| are taken into consideration. In addition, the k,± integrations in the equations 
of LCDAs and ,^-moments can be completed analytically when the Gaussian-type wave function 
is considered. After fixing the parameters that appear in the wave function, the curves and the 
corresponding decay constants of the LCDAs are plotted and calculated for the charmonium and 
bottomonium states. The first three non-vanishing ^-moments of the LCDAs are estimated and 
are consistent with those of other theoretical approaches. 
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I. INTRODUCTION 



Light-cone distribution amplitudes (LCDAs) of hadrons are key ingredients in the de- 
scription of various exclusive processes of quantum Chromodynamics (QCD), and their role 
can be analogous to those of parton distributions in inclusive processes. In terms of the 
Bethe-Salpeter wave functions ^{ui,ki±), LCDAs (f){ui) are defined by retaining the mo- 
mentum fractions Ui and integrating out the transverse momenta ki± p,]. They provide 
essential information on the non-perturbative structure of the hadron for QCD treatment 
of exclusive reactions. Specifically, the leading twist LCDAs describe the probability am- 
plitudes to find the hadron in a Fock state with the minimum number of constituents. In 
addition, the fact that 5-physics exclusive processes are under investigation in BaBar and 
Belle experiments, also urges the detailed study of hadronic LCDAs. In literature, there 
have been many non-perturbative approaches to estimate LCDAs, such as the QCD sum 
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, lattice calculation 
[qI, [lO|, Nabmbu-Jona-Lasinio model 



^chiral quark model from the instanton vacuum 
and the light-front quark model 



These studies have dealt with 
0, [l3, Q, axial vector Q, 0, 



1/CDAs of pseudoscalar 
15| . and tensor [6[ mesons. 
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The present paper is devoted to the study of leading twist LCDAs of p-wave heavy 
quarkonium states which include the scalar {xco^ Xbo), axial vector {xci, Xbi, hb), and ten- 
sor (Xc25 Xb2) mesons. The motivation of this study is as follows. Since the discoveries of 
J/ip and T, occurring more than thirty years ago, a great deal of information on heavy 
quarkonium levels and their transitions has been accumulated [l^. The numerous transi- 
tions between heavy quarkonium states are classified as strong and radiative decays, which 
shed light on aspects of QCD in both the perturbative and the non-perturbative regimes 
(for a recent review see [iTj]). In par ticular, some experimental results regarding Xcj mesons 



have recently been reported 



19|, 



20 



2l|. Therefore, a thorough understanding of their 



properties, such as LCDAs which are the universal non-perturbative objects, will be of great 
benefit when analyzing the hard exclusive processes with heavy quarkonium production. 

It is known that heavy quarkonium is relevant to non-relativistic treatments 22j. Al- 
though non-relativistic QCD (NRQCD) is a powerful theoretical tool for separating high- 
energy modes from low-energy contributions, in most cases the calculation of low-energy 
hadronic matrix elements has relied on model-dependent non-perturbative methods. In this 



2 



study, heavy quarkonium is explored within the hght-front quark model (LFQM) which is 
apromising analytic method for solving the non-perturbative problems of hadron physics 



231 as well as offers many insights into the internal structures of bound states. The basic 



ingredient in LFQM is the relativistic hadron wave function which generalizes distribution 
amplitudes by including transverse momentum distributions, and contains all the informa- 
tion of a hadron from its constituents. The hadronic quantities are represented by the 
overlap of wave functions and can be derived in principle. The light-front wave function is 
manifestly a Lorentz invariant, expressed in terms of internal momentum fraction variables 
which are independent of the total hadron momentum. Moreover, the fully relativistic treat- 
ment of quark spins and center-of-mass motion can be carried out using the so-called Melosh 
rotation [24]. This treatment has been successfully applied to calculate phe nomenologi cally 
many important meson decay constants and hadronic form factors 25|, l26|. 
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Therefore, the main purpose of this study is the calculation of the leading twist LCDAs of 
p-wave heavy quarkonium states within LFQM. 

The remainder of this paper is organized as follows. In Section II, the leading twist 
LCDAs of p-wave meson states are shown in cases of the vector and tensor currents. In 
Section III, the formulism of LFQM is reviewed briefly, then the leading twist LCDAs are 
extracted within LFQM. The ^-moments of these LCDAs are also calculated. In Section IV, 
numerical results and discussions are presented. Finally, the conclusions are given in Section 
V. 



II. LEADING TWIST LCDAS OF p-WAVE MESONS 



Amplitudes of hard processes involving p-wave mesons can be described by the matrix 
elements of gauge-invariant nonlocal operators, which are sandwiched between the vacuum 
and the meson states. 



(0|g(x)r[x, -x]qi-x)\HiP,e)), 



(2.1) 



where P is the meson momentum, e is the polarization vector or tensor (of course, e does 
not exist in the case of scalar meson), F is a generic notation for the Dirac matrix structure 
and the path-ordered gauge factor is: 

nl 



[x, y] = P exp igs ( dt{x - y)^A^{tx + (1 - t)y) 

" 



(2.2) 



This factor is equal to unity in the hght-cone gauge which is equivalent to the fixed-point 



gauge, {x — y) ^A^^ {x — y) = 0, as the quark- ant iquark pair is at the light-like separation |31 |. 
For simplicity, the gauge factor will not be shown below. 

The asymptotic expansion of exclusive amplitudes, in powers of large momentum transfer, 
is governed by the expanding amplitude Eq. fl2.ip . shown in powers of deviation from the 
light-cone x"^ = 0. There are two light-like vectors, p and z, which can be introduced by: 

p2 = 0, z'^ = 0, (2.3) 



such that p ^ P in the limit Mjj and z ^ x for = 0. From this it follows that {4] 

1 



Mi 
x'^ 



Px - J{Pxy - xmfj 



p 



M2 



_ (2.4) 



2Pz 



where Px = P ■ x and Pz = pz = y {PxY — x'^Mjj. In addition, if the meson is assumed 
that it moves in the positive 63 direction, then p"*" and z~ are the only nonzero component 
of p and 2, respectively, in an infinite momentum frame. For the axial vector meson, the 
polarization vector is decomposed into longitudinal and transverse projections as: 

e[ = -ip^-z^^\, e^^ = e^-el, (2.5) 
" pz \ 2pz J " 

respectively. For the tensor meson, the polarization tensor is: 

e^^im) = {ll]m'm"\ll]2m)e^{m')e''{m"), (2.6) 

(m is the magnetic quantum number) and e^,(= e^^ z^) can also be decomposed into longi- 
tudinal and transverse projections as: 

f" ( /If 2 \ 

= ^—[j,^^ - z^'^\ ^ e/;. _ _ (2.7) 

" pz \ 2pz J " 

LCDAs are defined in terms of matrix element of nonlocal operator in Eq. fl2.ll) . For the 
scalar (S), axial vector {A), and tensor (T) mesons, the leading twist LCDAs can be defined 
as: 

,M2 



{0\qiz)rqi-z)\SiP)) = fs / du e'^^^ 



pys{u) + z^-^gs{u) 
2pz 

, €Z 



(2.8) 



{0\q{z)Y'y5q{-z)\A{P,ex=o)) =ifAMA due'^'P^[p^—ct>A\\{u)+t^^gA±{u 



pz 

ez 



(0|g(^)7'^g(-^)|T(P,eA=o)) = /tM| I'du e'^^^ {p>^ ^ct>T\\{u) + '^qt^ 

Jo ^ [pz) pz 



:M^gT3{u)}, (2.11) 



2{pzf 



{0\q{z)a^''q{-z)\T{P,ex=±,)) = if^Mr / du e'^^^i^ 

+ (p-.^ -p^.-)^/.^||(«) + (6-.^ - 6r.-)^/.^3(^^,12) 

where u is the momentum fraction and ^ = (1 — m) — « = 1 — 2u. Here 05, 0a,t||5 and 
4>A,T± are the leading twist-2 LCDAs, and the others contain contributions from higher- 
twist operators. Due to G-parity, 0s, gs, 03Ai±, hsA^w, hsA^s, 0iai||, g^A^±, g^A^s, 0t||, gT±, 
gT3, 4>T±, hxw, and Hts are antisymmetric (odd) under replacement u 1 — u, whereas, 
0iyii±; ^Miih ^iyii3; 03Ai||; g^AiL, ^^d (y's^jS are Symmetric (even) in the quarkonium states. 
Therefore, the leading twist LCDAs are normalized as: 

1 p1 

(odd)/ \ _ ^ / J„,J,(even), 



du^(j)^°'"'>{u) = 1, du(f)^'''"''>{u) = 1. (2.13) 

Jo 

and can be expanded ^ in Gegenbauer polynomials C^^^(0 as 



</'(^,/^) = 0as(O 



Li=0 



(2.14) 



where 0as(O = 3(1 — ^^)/4 is the asymptotic quark distribution amplitude and ai{fi) are 
the Gegenbauer moments which describe to what degree the quark distribution amplitude 
deviates from the asymptotic one. Cf^^(^)s have the orthogonality integrals 

£(1 - e)criocTiOd^ = (2-15) 

Then a/ can be obtained by using the above orthogonality integrals as 

An alternative approach to parameterize quark distribution amplitude is to calculate the 
so-called ^-moments 

{n, = /_/er0(e,/i). (2.17) 



To disentangle the twist-2 LCDAs from higher twist in Eqs. fl2l8D ~ fl2TT2l) . the twist-2 
contribution of the relevant nonlocal operator q{z)rq{—z) must be derived. For the F = 
7'^ (75) case, the leading twist-2 contribution contains contributions of the operators which 



are fully symmetric in Lorentz indices 



32 



33|: 



[<l{-zm^M^)h = E -^I0)p-f + ^ D^AilMO), (2.18) 

^ n! I n + 1 n + 1 J 

where D = D — D and D = d — igB"'{X"' /2). The sum can be repressed in terms of a 
nonlocal operator, 

[qi-z)YhMz)h = C dt^q{-tz) AlMtz). (2.19) 

Jo oz^ 

Taking the matrix element between the vacuum and the p-wave meson state, we obtain: 

{0\[q{-z)Yq{z)]2\SiP)) = fs ^/w05(«)|p'^e^^^^ + (P'^ - p'^) cite^^*^^ J(^.20) 



{Om-z)rjrMz)],\A{P,ex=o)) = iIaMa /_ du<PA\\{u)\p^-e'^^-' 

I pz 



+ -P^^) ^^c/te^«*P"|, (2.21) 



(0|[g(-z)7^g(^)]2|T(P,eA=o)) = frM^ du<PT\\{u)L 



p__ 



{pzy 



/ urn mm \ „i 

+ 2 P^^^ / dte'^'P' ;2.22) 

\ »2; (vzV Jo 



pz (pzy J Jo 

The derivations of Eqs. (12.201) ~ (12.221) are shown in Appendix A. We can use Eq. (12.181) . 
and then expand the right-hand sides of Eqs. (12.201) ~ (I2.22p . as 

^ 1 ((z-Dr .. niz-DY-^ 



^ n! I n -h 1 n -f 1 J 

fsT.^ du4>s{u){^pzr\p'' + iP'-Pn dte\, (2.23) 

n=0 I JO ) 

E A(0k-(0){^^1^7'^ + ^^^^:^D'^4g(0)|A(P,6)) 
ln + 1 n + 1 J 

'/.M,|^/;<;..^,„(.)Kp.)"{p''|+ (2.24) 

f^M^J2- du^T\\{u){^pzr{p''--^ + 2i P'r^] dtr\, (2.25) 

^n!Jo [ [pzy \pz [pzy J Jo J 
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respectively. Picking n = in Eqs. (12.231) and (12.241) . we obtain 



(0|g(0)7^g(0)|5(P)) = fsP^ / du4>s{u), 

Jo 



(2.26) 
(2.27) 



Note that the tensor meson cannot be produced by the V — A current, then we pick n = 1 
in Eq. (I2.25P and obtain 

^(0|g(0)(7^^-:D+ /D'^)g(0)|T(P,eA=o)) = frM^e^^' f' du^4>Ti\{u). (2.28) 
2 Jo 

From the normahzation Eq. f|2.13D . we have (0|g7'^75gpy4i(P, e)) = ifzA^MiA^^^ which is 
consistent with the results of |34| . 

Next, we consider the case of F = (7^,^(75), where the leading twist-2 contribution contains 
contributions of the operators: 



n=0 



2n + l 2n + l 



2n + l 

The sum can be also represented in terms of nonlocal operators: 



^V^" (75)g(0). 



(2.29) 



Jo 



dz, 



'dz^ 



Taking the matrix element between the vacuum and the axial- vector and tensor meson state, 
we obtain: 



{Om-z)a^''j,q{z)h\A{P,e,=^,)) 









(j)A±{u) 







{Om-z)a'^'^qiz)],\TiP,e,=^,)) 



pz 



(2.31) 



pz 



pz 



+ {hT\\{u) - 0T±(ti)) T'^"e'«P" + (^^^ - 3T'^'') j^^ dte'^''^' |, (2.32) 



where 



(e^P" - e'^P^) - {e^z" - e^z 



2pz 



1 



/Vf2 

pz 



2pz 



pz 



(2.33) 



The derivations of Eqs. (12.311) and (I2.32p are shown in Appendix B. In contrast to Eqs. 
(I2.20p ~ (12.221) . the twist-2 LCDAs do not disentangle entirely from the higher twists in 
Eqs. (I2.3ip and (12.321) . Taking the product with e^p,Zy and e^p,,Zy in Eqs. (I2.3ip and (12.321) . 
respectively, to obtain, 



-^)^^'^'e±/.75g(^)]2|A(P,eA=±i)) = /i t du<^Ai.{n)ht-t^Pz) 

Jo 2 



{0m-z)af''e^„qiz)]2\T{P,ex=±i)) = if^Mr / du(t)T±{u)-e^'e^„ 



Jo 

(2.34) 

Jo 

(2.35) 

Then, we use Eq. (12.291) and expand the right-hand sides of Eqs. (I2.34p and (I2.35P as: 

£1(01^(0) ^^VJ^';^^" ^"'^^M75g(0)|A(P,6,=^0) 



n=0 



2n + 1 



oo -n ^1 1 

^0 ^' -^0 2 



1+ / dtt 





2n 



(2.36) 



E :^(0|glO)^^^±il^^^a^'-6^,g(0)|T(P, 6.=±0) 



n=0 



2n + l 



oo -n ^1 

E' 

n=0 

Picking n = in Eqs. (I2.36P and (12.370 . we obtain: 



oo -n ^1 1 

if^MTY.-J du<PT±{u)-e^'e^,.{^pzy 
nl Jo 2 



1 + / dtt 





2n 



(0|g(0)a'^'ex^75g(0)|A(P,eA=±i)) = fi f' du(f)A±{u){e ■ e^Pz), 

Jo 

(0|g(0)a^*e^^.g(0)|T(P,eA=±i)) = tf^Mr f du(l)T±{u)e^'e±,.. 

^ 



It is worth noting that the author of Ref. 



(2.37) 



(2.38) 
(2.39) 



3l| also considered an approach that disentangled 



the twist-2 LCDAs from the higher twists, in the case of an axial vector meson state (F = 
(7^*^75): Besides Zy, Eq. (I2.3ip has taken the product with a term proportional to {e^Pz — 
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Pfj,ez). We find this approach equivalent to ours. The derivation is as follows. The term 
{e^j^Pz — Pfj,ez) can be expanded by using Eqs. (12. 4p and (12.51) as 



ez „ ezMl ezM\ 



^uPz — Pu^Z = Pa Pz — Zu—, ZT^Pz + e^nPz — Pu^Z — Zn 



ezM' 



-z 



pz 



+ e^^Pz. (2.40) 



The first term of last line has no contribution to the result because a^'^ is antisymmetric. 



III. GENERAL FORMULISM IN LFQM 



A. Framework 



A meson bound state, consisting of a quark qi and an antiquark q2 with total momentum 



P and spin J, can be written as (see, for example 26|): 



|M(P,2^+i Lj, J,)) = J {rf3fci}{t/3fc2} 2{2nfS\P - ~h - h) 



X E ^i^1^i,^2,Ai,A2) |gi(A;i,Ai)g2(^2,A2)), (3.1) 



where ki and /c2 are the on-mass-shell light-front momenta. 



and 



~k={k+,k^), k^ = {k\e), k- = '^'^^ ^ , (3.2) 



{d-^k} = 



2(27r)3 ' 

\q{k,,X,)q{k„\2))=blik^)dl{k2m, (3.3) 

{by{k'),b{{k)} = {dy{k'),d{{k)} = 2(27r)3 6''Ck' - k) 6yx. 

In terms of the light-front relative momentum variables {u,k±) defined by 

fc+ = (1 - m)P+, fc+ = mP+, 

ki± = (1 - u)P± + K±, k2± = uP± - K±. (3.4) 

The relative momentum in z direction can be written as 

_uMo ml + Kl 
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The momentum-space wave-function ^ls' ^ Lj meson can be expressed as 

■^ifikuh,Xi,X2) = -^{LS;L,S,\LS;JJ,)RfXiu.^^±) ^llA%^^±), (3-6) 



where 'PllA'^^ '^i) describes the momentum distribution of the constituent quarks in the 
bound state with the orbital angular momentum L, {LS; L^SzlLS; JJ^) is the corresponding 
Clebsch-Gordan coefficient and R\^x^ constructs a state of definite spin (5, Sz) out of light- 
front helicity (Ai,A2) eigenstates. Explicitly, 



where |sj) are the usual Pauli spinors, and TZm is the Melosh transformation operator 25l |: 



[s\nM[Ui,K^,mi)\\) = ; — , (3.8) 

{rrii + UiMoY + 4 



with ui = 1 — u, U2 = u, and n = (0, 0, 1) is a unit vector in the ^-direction. In addition, 

= (ei + = + HilA^ (3.9) 

Ui U2 



where Mq is the invariant mass of qq and generally different from the mass M of meson 
which satisfies = P^. This is due to the fact that the meson, quark and anti-quark 
cannot be simultaneously on-shell. We normalize the meson state as 

(M(P', J', J'z)\M{P, J, Jz)) = 2{2nf P+6%P' - P)6rj5j,j^ , (3.10) 

in order that: 

^^^^y 'Pl'l'S'^^ i^a.)vllAu, = h',L h'^,L.- (3.11) 

Explicitly, we have: 

V^iL. = f^L.^p, (3.12) 

where kl^=±i = Tii^±x Hf^Ly) / V^, kl^=o = f^z are proportional to the spherical harmonics 
YiL^ in momentum space, ifp is the distribution amplitude of p-wave meson. In general, for 
any function F{k,), (pp{u) has the form of: 



Mu) = nJ^F{k), (3.13) 



10 



where the normahzation factor is determined from Eq. (13.110 . 



In the case o 



tdSSz 



25|,l28|,l36|: 



a p-wave meson state, it is more convenient to use the covariant form of 



"-1 "-2 



^2 Mo(Mo + mi +m2) 

XM(fci, Ai)(^ + MQ)V2s+ip,v{k2, As), (3.14) 



where 



Mo = - (mi - m2)2, 

u{k,X)u{k, A') = 

2777, 

v{k,X)v{k, A') = -— 5a,a' 



P = fci + ^2, 
u{k, \)u{k, A) = 

A 

J2vik,X)v{kA) 

A 



^ — m 
^ k+ ■ 



For the scalar, axial-vector, and tensor mesons, we have: 

If K-P\ 



5; 



rap, = e,,ri-Kn, 



e-i^ 75, 



(3.15) 



(3.16) 



where K = {k2 — ki)/2 and: 



tA=±l 



— e-l(±l) ■ P^, 0, e-l(±l) 



el(±l) = T(l,±2)/V2 



tA=0 



1 (-M^ + Pl 



(3.17) 



Mo V ^+ 

Note that the polarization tensor of a tensor meson satisfies the relations: e^i, = e^^ and 
t^yP^ = = 0. Eqs. (I3.14p and (13.161) can be further reduced by the applications of 
equations of motion on spinors: 



{lS;L,S,\lS;JJ,)kLzRfX{u,f^±) = ^^^j^uik,, Xi)ns+ip^v{k2, X2), (3.18) 

v2 Mq 



where 



2V3M0' 
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r'3 



3Pl 



Tip 



e • ^75, 

'Ij- ^Ml - 2e ■ K{7m - ma)) 75, 
2K^' 



2V2M0 



Mo + mi + m2 



(3.19) 



B. Analysis of Leading twist LCDAs 



Next, the matrix elements of Eqs. dOi), dMZD, (E3SD, and fl^:^ will be cal- 

culated within LFQM, and the relevant leading twist LCDAs are extracted. For the scalar 
meson state, we substitute Eqs. (13.11) . (13.61) . and (I3.18P into Eq. (I2.26P to obtain: 



m2Yqi\S{P)) = Kj{d'k,} ^ifih,k2,Xi,X2){0\q2Yqi\qiq2) 

fi + mi \ Mq 



Al,A2 



r 



/(fa + ma' 



fsP^J duct>{u). 



(3.20) 



For the "good" component, /i = +, the leading twist LCDA 05 can be extracted as: 

v^2iVc r cPkx. [(1 - M)m2 - Mmi]Mo 



05(m) 



fs J 2(27r)3 ^3m(1 - u)Mo 
In the quarkonium case mi = ma = m, Eq. (13.211) can be further reduced as: 

V2 f d^K± (1 - 2M)m 
(ps{u) = -J- / -7-—^^===iPp{u,K^). 
Is J 2{2ny Ju(l-u) 



(3.21) 



(3.22) 



A similar process can be used for the axial vector and tensor mesons which correspond 
to Eqs. (E27D, fl238l) and (ESHD, flCTl) . respectively, and the leading twist LCDAs are 
extracted as 



03Ai|](u) 
0iAi||(m) 



2^3 f d^K± 



f.A, J 2(27r)3 _ u)Mq 

a/6 r d'^K± (1 — 2'u)m 



f^A, J 2(27r)3 y^IT^ 
ye r d'^K^ (l-2n) 



03Ai±(m) 



h J 2(27r)3 

73 r d'^K^ (1 - 2M)m 



/4 ^ 2(27r)3 y^IflT^ 



V9p(M, K_l) 



Mn-m- 



Mo + 2m 



(3.23) 
(3.24) 

^Pp{u,K±), (3.25) 
(3.26) 
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0iAi±H = 7^ / oro^^s r— — <^p(m,«:±), (3.27) 



V6 /• d^ft;^ (1 _ 2 



u] 



J 2(27r)3 ^ 



m + 



2k' 



2 



Mo + 2m 



K^). (3.28) 



From the normalization Eq. ( ]2.13p . some relations among the constants /a/s could be easily 
obtained as: 

V^/S = /Mi = V^/sii = /odd, and ^ = /i^^ = /even, (3.29) 

then, the relations among the relevant LCDAs are: 

0s = 0iAi|| = 03Ai± = 0odd, and 03Ai|| = = 0even, (3.30) 

where the subscript "odd (even)" means the odd (even) function of u and 

0odd(M) = / o/o ^3 r—^ =VlAU,Ki_), (3.31) 

/odd 2(27r)3 ^'^(i _ u) 

0even(M) = / :T77^T5^=^= </'?(«, /t±) (3.32) 

/even 2(27r)'^ 1 - M Mn 



Note that Eqs. fl3.29p and (13.301) are independent of the form of ipp. Furthermore, in the 
nonrelativistic situation, the momenta Kz^± are much smaller than quark mass m, and Mq 
can be reduced to approximately 2m. Thus, we have: 

/odd ^ /t ^ /t , 0odd ^ 0T|| ^ 0T±- (3.33) 

and 

... V6 /■ c?V (l-2?z) Mo 

0odd(M) ^ / — — ^^== — V2p(M,«± . (3.34) 

/odd 2(27r)3 1 _ w 2 



From Eqs. (13.321) and (13.341) . one can obtain: 

/odd — /even, (3.35) 

and relate the ^-moments of the 0even to those of the 0odd as 



13 



with the function F = F{\k,\). The derivations of Eqs. 
Appendix C. The above results are consistent with those of 



[3.351) and (13.361) are shown in 



6| in the nonrelativistic approx- 



imation. 

Next, we choose a Gaussian- hke wave function, as shown in 



28|: 



4^2 / TT \3/4 dK, 



\K\ 



(3.37) 



for further calculations. In Eqs. (I3.3ip and (13.321) . the k± integrations can be performed as 
follows: 



4>odd{u) 
4>even{u) 



^/6(l -2u)m(2 



/odd 
/even \ ^ / 





5/4 


[5 1 












4'" 





r 1 1 


+ 3r 


[3 1 








1 








4'^. 





(3.38) 
(3.39) 



where w = d/[4:u{l — u)], d = m^/2/3^ and: 

r[a, w] = 



is the incomplete Gamma function. The incomplete gamma function may be expressed quite 
elegantly in terms of the confluent hypergeometric function: 



r[a, w] = T[a] — a w°' x iFi(a; a + 1; —w) 



(3.40) 



where 



(3.41) 



^0 (ai)„(a'2)„...(a^)„ n\ 

and {a)n = {a + n — l)\/{a — 1)\ is the Pochhammer symbol. However, the k± integrals in 
Eqs. (I3.25P and (I3.28P cannot be analytically performed. The crux is the term proportional 
to 1/Mq + 2m. This term may be rewritten and expanded as: 



Mo + 2m 



4m (l + 



Mo-2m 
4m 



4m 



Mo -2m\ 



4m 



(3.42) 



^ For the tensor mesons, Ref. 16| has the relation fx = y |/t- The additional y | factor is the Clebsch- 
Gordan coefficient of the polarization tensor e^^ for the tensor meson state T{P, e\=o)- This distinction 
is from the different definition of ^||: There is a e^i, in both hand sides of Eq. (|2.1ip . By contrast, in 
the upper part of Eq. (4) of Ref. [6|, e^i/ appears only in the left hand side. 
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We only consider the first two terms in the square bracket of Eq. ( 13.420 because Mq— 2m goes 
to zero in the non-relativistic approximation. Then, the approximate form of, for example, 
(w) which containing the first one term and first two terms in the square bracket of Eq. 
f l3.42p are defined as: 



V6 f (Pk^ (1 - 2u) 



It J 2(27r)3 ^^1(13^ [ 
r cPk^ (1-2m) 



Mo-m- 

4m 



Mo - m 



fr J 2(27r)3 _ ^) 
respectively. The k± integrals can be performed as: 
VQ{l-2u)mf2 



4m 



Mo - 2m 
4m 



(3.43) 
V?p(m, k^),(3.44) 



/t 

VQ{1 - 2u)m f 2' 









[5 1 




[9 1 




)""\ 


















4w 










[5 1 


-±T 


r9 ] 




















4'^. 


8w 


4'^. 



+ ^=^r 

v/2m 
15/5 , 



[7 1 






1 







(3.45) 



+ 



4^2; 



m 



7 ■ 
4'^. 



11 



(3.46) 



4:\/2wm 

We can find that the both curves of 0^|j(-u) and 0y|j(M) are almost overlap each other by 
using the parameters in Table 1 (See Sec. IV). Therefore we only consider the first term in 
the square bracket of Eq. fl3.42p and take the approximation ^ 0y||(M). Finally the 

form of (f)T±{u) is: 



V6{1 - 2u)m (2 



5/4 



2/ 



e'^ r 



T 



7r 



4' 



W 



,W 



(3.47) 



In addition, the ^-moments of Eqs. fl3.38p . fl3.39p . (13.451) . and (13.470 can be analytically 
expressed as: 



A. 



oddS - 





-5" 
4 


r 







"i" 

4 


r 





r 


5' 
4 




r 


1 + / 



/: : —d 

2 ' 4' , 



(3.48) 



2-^2 



3 _^ 3__1 7__^ 

"4 " 



3r 


"-1' 
_ 4 




4r " 







-5" 
4 _ 


r 


4 ^ f 





3 -7 - ^; 7; + 2i"2 -7 - 



(3.49) 



9 5 13 
4'4'T' 



15 



+ 



V2T 


4 


r 





|_.f.(_|_,4_,),(9,40.f.(-f-<;-i;-.; 
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P 



m 



r 


'7' 
4 




r 


l + i 



3 3 \ d'/'T 



+ 1 



(3.50) 



A 



4 ^ 



^ , 1 ,9 A 5^/1 , 9 5 13 



2r 

1 


"9 
4 




dT 


l + l 



where 



^odd 

At 





5/4 


[3 ,1 




) e^r 











A^ 





5/4 


3 , 




) e'^r 


- + / 






2 



A 



T± 





5/4 






(tt; 


























^5/4 















and / is a non-negative integer. The derivations of Eqs. (13.481) ~ (13.511) used the formula: 



iFi(a;6; c) 



6-1 



[iFi(a;6-l;c)- iFi(a - 1; 6 - 1; c)] 



(3.52) 



which is easily checked from the definition of the confluent hypergeometric function Eq. 
IV. NUMERICAL RESULTS AND DISCUSSIONS 



In this section, the LCDs, constants /s, and ^-moments are estimated. Prior to numer- 
ical calculations, the parameters m and (3, which appeared in the wave function, must be 
determined firstly. We consider the Hamiltonian of the p-wave heavy quarkonium state as: 



H = 2V m? + K? + hr 



3r 



+ giS ■ L + g2Si2 + gzSi ■ S2, 



(4.1) 



where br (— ^) is the linear (Coulomb) potential, S12 = (3si ■ f S2 ■ r — Si ■ S2) is the tensor 
force operator, and (71,2, 3 are the functions of the relevant interquark potentials (the details 



are shown in, for example. 



37 



38|). In this way, the mass difference between the spin- 
single ground state M(^Pi) and the spin-weighted average of the triplet states Mi^Pj) = 
[M(3Po) + 3M(3Pi) + 5M(3p2)]/9 only has the contribution which comes from the spin-spin 
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interaction.^ Experimentally this hyperfine splitting is less than 1 Mev [16| in charmonium 
sector, and can be neglected here. Then, we can use the mass M{^Pj) and its variational 
principle for the Hamiltonian Eq. fl4.ip in order to determine parameters m and f3. In the 
process, the conjugate coordinate wave function of Eq. fl3.37p 



3 7rV4 



/5r exp 



(4.2) 



is necessary. The values of the additional parameters b and come from literature [39| : 



6 = 0.18 GeV^ 



0.36, 



(4.3) 



for the heavy quarkonium states. We individually vary b and to realize how they connect 
to m and (3. The results are shown in Table I. We find the parameters m and j3 insensitively 
depend on b and ag. 



TABLE I: The connections among the parameters b, as and m, f3 of the p-wave heavy quarkonium 
states. 



6(GeV2) 


as 


mc(GeV) 


/3cc(GeV) 


m;,(GeV) 




0.18 ±0.02 


0.36 


1 00-0.03 
J^-J"+0.04 


0.4891°:°}! 


4.76 =F 0.02 


r, 701+0.022 
U-'y-L-0.023 


0.18 


0.36 ± 0.04 


1.38 ±0.01 


0.489 ± 0.007 


4.76 ± 0.02 


r, 701+0.022 



Next, we use the center values of the parameters in Table I to calculate and plot the 
LCDAs in Eqs. (E3HD, dDD, (IS3SD, and (Km . The results are as follows: 

/odd = 0.0884 GeV, 
/t|| = 0.124 GeV, 

for the charmonium states, and 

/odd = 0.0674 GeV, 
/y|l = 0.0750 GeV 

for the bottomonium state. The curves of LCDAs are shown in Figs. 1 and 2. These 

^ The calculations of the expectation values of the fourth and fifth terms for the Hamiltonian Eq. (14. 1|) can 
refer to the appendix of Ref. [4^] . 



/even = 0.109 GcV, 

/^^ = 0.0978 GeV, (4.4) 

/even = 0.0716 GcV, 

/^^ = 0.0692 GeV, (4.5) 
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FIG. 2: The LCDAs 6 

cven(0 of the charmonium (sohd hne) and bottomonium (dashed hne) states. 

results are consistent with Eqs. (13.331) and (I3.35p . In addition, the curve of 4>^^{^) in Fig. 
2, in which ^ is peaked around zero, was sharper than that of (f)'^'^{C,). This meant that 
the momentum fraction u in the bottomonium state is more centered on 1/2 than in the 
charmonium state, which is reasonable, as the mass of b quark is larger than that of c quark. 
A similar situation exists for the odd functions in Fig. 1. 

Finally we show the LCDAs in terms of the ^-moments. Eqs. (I3.48p . (I3.50p . (I3.5ip . and 
(I3.49P are calculated for Z = 3, 5, 7 and I = 2,4, 6, respectively. The results, which compare 
with the other theoretical evaluations, are as shown in Tables 2 and 3. In Table 2, Refs. 
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6l used the QCD sum rules with the non-relativistic wave functions. The authors of Ref. 

4l| calculated in the framework of the BuchmuUer-Tye potential model, and found their 
results are in agreement with experiments which included the leptonic widths and hyperfine 
splittings. The authors of [42] calculated in the Cornell pot ential. For the charmonium 
sector, our results are not only consistent with those of g], 41; 42], but also conform to Eq. 



fl3.36p . In the framework of NRQCD, the authors of Ref. j6|] related the relative velocity 



TABLE II: The ,^-moments for the p-wave charmonium states. (C*)0odd ~ iOhc- ^ iO 



+ 1)) 



moTTipnt, 


this work" 


[61 tt 


[411 


[421 


/<Podd 


n 1 QD 


n 1 a _|_ n no 


n 1 8 


1 6 


('?^)</'odd 


0.0507 


0.050 ±0.010 


0.047 


0.040 


('?'^)<^odd 


0.0164 


0.017 ±0.004 


0.016 


0.013 




0.206 










0.0583 










0.0198 










0.188 










0.0498 










0.0160 








)<^cvcn 


0.0662 


0.06 ± 0.01 






)<^cvcn 


0.0110 


0.010 ±0.002 






)(/)gvcii 


0.00261 


0.0024 ± 0.0006 







TABLE III: The ^-moments for the j?-wave bottomonium states in this work. 



)(/'odd 


0.0666 




0.0691 




0.0665 


)0cvcn 


0.0226 


(^^)0odd 


0.00685 




0.00735 




0.00684 




0.00142 


(?'^)0odd 


0.000922 




0.00102 




0.000919 




0.000139 



of quark-antiquark pair inside the p-wave charmonium state to the ^-moments as: 

« -I- 3 

{v-), = ^{C^') + Oiv-+') (4.6) 
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If the ^-moments (^")</,odd considered, we obtain: 




0.317, 



0.118 



0.0492, 



for the charmonium sector and 




0.111, 



0.0160, 



0.00277, 



for the bottomonium sector. These results are consistent with the values (f^)cc ~ 0.3 and 
(f2)feb ^ 0.1 used in NRQCD. 

V. CONCLUSIONS 

This study discussed the leading twist LCD As of the p-wave heavy quarkonium states 
within the light-front approach. The twist-2 LCDAs have been disentangled from the higher 
twists by appropriately coping with the nonlocal operators q{—z)Tq{z). For the T = 0-^,^(75) 
case, we proved that our method is equivalent to that of Ref. [3l|. Next, these LCDAs have 
been shown in terms of the light-front variables {u, k,±) and the relevant decay constants. 
We found that the decay constants and LCDAs had the following relations: -\/3/s = f^A^ = 

^fUii= /odd), fujV^ = /iii(= /even), and 0s = = 03^^_l(= 0odd), (p^A^W = 4>iAi±{ = 

0even)- If oue takcs the non-relativistic limit and the wave function as a function of \k\, 
then the above relations among the decay constants could be further simplified as /odd — 
fx — fx — /even, and in addition, the ^-moments of 0odd and 0even have the relation: 



The integrations for the equations of LCDAs and ,^-moments could be analytically 
performed when the Gaussian-type wave function is considered. The parameters m and /3, 
which appear in the wave function, were determined by taking the mass of the spin-weighted 
average of the triplet state Mi^Pj) and the variational principle for its Hamiltonian into 
account. We found the parameters m and j3 insensitively depended on the linear potential 
constant h and the strong coupling constant a^. The curves and the corresponding decay 



(n<^o../(^ + i)- 
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constants of the LCDAs 0odd, (pT\\, 0t±, and 0even were plotted and calculated for the 
charmonium and bottomonium states. These results are consistent with the relations which 
are mentioned^in last paragraph. However, the value of /odd is about a factor of two smaller 
than that in [6] which was studied within QCD sum rules. In addition, the first three ^- 
moments were calculated, and were consistent with those of other theoretical approaches. 
The relative velocity of quark- antiquark pair (f^) of charmonium and bottomonium states 
were also estimated and were consistent with those used in NRQCD. 
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APPENDIX A: DERIVATIONS OF EQS. (f2:20]> - ([2722]) 



Firstly, Eqs. and fl2llD are rewritten as 

{0\q{z)j,q{~z)\S{P)) = fsj^ du e'^^UP.^siu) + z.^Jgsiu) - 4>s{u)n , (Al) 



{0\q{zhf.l5q{-z)\A{P,ex=o)) = iJaMa j du e'^^'{e^(t)A\\{u) + 



tz 



'2{pz) 



e±M [9a±{u) - (j)A\\{u) 
gA3{u) - (j)A\\{u)] }(A2) 



\q{zh,qi-z)\TiP,ex=o)) = ItM^t / rfw e^«P^{^0T||(«) + ^ [^7r±H - <^t||H 



}(A3) 



respectively. Next, we sandwich both sides of Eq. (12.191) between the vacuum and, for 
example, the scalar meson state 

{m{-^)iMms{p)) = I' dt-^m-tz) Mt^)\siP)) 

1 /~) 1 

= fs [ dt—Pz [ due'^^'Mu) 
Jo ozi^ Jo 

= fs du(t)s{u){p, dte'^'^P^ +p^{iiPz) jytte'^'A.{M) 



The second term of the last line can be further calculated as: 



pMPz) / dtte'^"^' =pa — / dtt 
Jo pz Jo 



rz I . . . a_ .^ip^ _ 



Jo 



(A5) 
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We can substitute Eq. (lA4p for Eq. (lASP and obtain Eq. fl2.20p . In addition, the same 
process can be used to obtain Eqs. (12.211) and (12.221) . In fact, the above process has been 
used for the vector meson state in Ref. |32j. 



APPENDIX B: DERIVATIONS OF EQS. ([231]) AND ([232]) 



Eqs. (12.1 op and (12.120 can be rewritten as: 



J 



u) 



Ml 



(0|g(z)a^,g(-^)|T(P,e,=±0) = ^/t^t l'due'^P'{ 



(Bl) 



^IJ,»Pl/ ^•uPfJ.) , / N 

0T± [U) 



V^,[hT.{u) - <PTdu)]l (B2) 



respectively. Then, we sandwich both sides of Eq. (12.300 between the vacuum and, for 
example, the axial-vector meson state, 



(0|[g(-^)7M75g(^)]2|A(P,eA=±i)) 
1 



dt 



±-{Q\q(-ez)a,,^,q{ez)\A{P,,))+z-^^^^ 



/a / du\(t)AL{u) 



ie,P,-e,P,) j\te'^'"P' + 2pzS,M) dtt'e'^'"^' 



+ {hA\\{u) - (t)AL{u)) 



, dte'^'-P' + 2pz%M) / dtt'e 
Jo 



(B3) 



We can further calculate the integral as: 

1 /-i (9 



e / dtt^e'^^'P' = — I dtt^e'^''^' = — e'^P' - f dte'^''^' 
Jo 2pz Jo at 2pz I Jo 



(B4) 



2pz Jo at 2pz 

and then substitute Eq. ( IBSp for Eq. ( ]B4p to obtain Eq. (12.310 . The same process can be 
used to obtain Eq. (I2.32p . 
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APPENDIX C: DERIVATIONS OF EQS. (l335]> AND (l336]> 



From Eqs. (I3.32p . (13.341) and the normalization Eq. fl2.13p . we have 

£k± ^ Mo 



/odd ^ V6 1 '^rf^ / 

/even = j j 



2(27r)3yr^F 2 
(i^K^ 1 



(CI) 
(C2) 



2(27r)3v/r^Mo 

Taking Eqs. (ICip and (lC2p integration by parts with respect to n\ and ^, respectively, and 
using normalization conditions: v^p(^, k±) must go to zero when ^ (k_l) go to ±1 (infinity), 
and thus, we can obtain: 



/odd ^ -VeJ'^d^^ J 

/even = J^^d^^ J 



d'^Kx. 



i d 



2(27r)3 



Mo 



(C3) 
(C4) 



From Eq. (13.130 . the differentiations in Eqs. (1C3P and (IC4p can be expanded as: 



i d 

VT^dKl 
d_ 
di 



Mn 



v/r^Mo 



f(k) + v.. . . 



F(/«) + 



2(i-e)c?4 

1 



(C5) 
F(/?), (C6) 



2VMo(l - ^2)2 V ^ ^(1 _ ^2) rf^ 

respectively. If the function F = F{\k,\), one can expand the arbitrary function -F(|k|) as a 
polynomial of \k\: 



Fi\K\)=Y.cM'- 

s 



(C7) 



Using the relation \k\ = JM^/A — m?, one can calculate the differentiations in Eqs. (ICSP 



and (C6) as: 



iMl'^ d 



^(W) = 7SaT^^^(W) = ?=.i(r3e)i(^-'"j ■(«) 



2(1 -e)tiK2 " ^' " VM, 
Then, Eq. (I3.35P is obtained. Regarding the ^-moments of 0odd and ^even, we have: 



)</'od 



v/6 

/odd -^-1 

/ever 



n+l 



rf2/t_L e Mo 



2(27r)3v/r^F 2 



d^K^ 



2(27r)3^/r^Mc 



(C9) 
(CIO) 
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Using the above processes, Eq. fl3.36p can be obtained. 
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